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Ionizing  radiation  studies  involving  210  rhesus  (Macaca  mulatta)  monkeys  has  been 
t'xnr.ined.  These  studies  contain  n wide  variety  of  experimental  parameters  pre- 
l art'd  ly  investigators  interested  in  unique  situations.  When  proposed,  they  were 
viewed  as  pilot  studies  to  determine  the  presence  or  the  absence  of  radiation 
effects  upon  performance.  This  report  examines  the  studies  a posteriori  grouping 
them  as  follows  : Distribution  Gpecif  Ication  , ;EDrl(^  Determinations , Dose  Rates, 
Mixed  Rates,  and  Descriptive. 
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Distribution  specif icution  studies  found  that  the  onset  time  of  each  episode 
had  an  exponential  distribution,  and  the  duration  of  each  episode  was  either  an 
exponential  or  Weibull  distribution.  The  presence  cf  similar  distributions 
across  many  experiments  suggests  that  the  same  median  lam  is  at  work  in  each  study 
and  episode,  differing  only  in  intensity.  LDr.p  studies  found  the  effective  lose 
for  SOI  emetic  response  to  be  *'» U 6 rads  (Co1’0,  CO  ruds/niin),  and  that  motion 
reduced  the  EDcjq  to  2r>P  rads.  Dose  rate  studies  indicate  a requirement  for  nt^ 
least  10  animals  per  treatment  group  in  making  intergroup  comparisons.  Mixed 
rate  and  descriptive  studies  indicate  the  need  for  more  systematic  investigations 
In  an  average  rhesus  monkey,  if  emesis  occurred,  the  first  episode  occurred  at 
^0  min  postirradiation,  followed  in  about  15  min  by  a second  episode,  and  in 
another  15  min  by  a third  and  last  episode. 
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INTRODUCTION 

This  report  summarises  the  contents  of  the  USAF  School  of  Aerospace 
Medicine  Radiation  Sciences  Division  ionizing  emetic  repository  from 
1971  to  1977.  This  period  includes  the  21  studies  listed  in  Table  1. 

As  a group,  the  studies  contain  a wide  variety  of  experimental 
parameters  prepared  by  different  investigators  interested  in  unique  and 
sometimes  specified  situations,  and  emesis  data  was  collected  incidental 
to  the  primary  objectives  of  their  experiment.  Nevertheless,  the  emesis 
data  can  be  viewed  as  pilot  studies  reflecting  on  this  common  syndrome, 
recognizing  that  the  experiments  were  not  designed  for  elaborate 
modeling  efforts.  With  this  in  mind,  we  grouped  the  studies  a posteriori 
on  the  basis  of  the  analysis  permitted  by  the  data  collected  as  follows: 

Distribution  specification  (Studies  2,  18,  19,  27,  30,  33) 

ED50  determinations  (Studies  39,  **l) 

Dose  rates  (studies  31.1,  31.2) 

Mixed  rates  (studies  36,  37,  38,  Uo) 

Descriptive  (Studies  3.1,  3.2,  26,  28,  32,  3^,  35) 

Computer  programs  used  to  examine  this  repository  are  documented  to 
provide  continuity  to  this  data  base.  Brief  explanations  are  given  so 
nonprogrammers  can  access  the  data  to  update  and/or  pursue  additional 

analysis . 

We  use  the  operational  definition  of  emesis  given  by  Fleischman  and 
associates  (k).  That  is,  an  emetic  episode  is  defined  as  a "series  of 
thre.  or  more  reflexive  movements  consisting  of  coordinated  mouth 
openings,  opisthotonus,  and  maximal  abdominal  musculature  contraction." 
Tnterepisode  time  is  defined  as  "3  or  more  minutes  between  the  last 
emetic  movement  of  an  episode  and  the  first  emetic  movement  of  the 
subsequent  episode."  Note  that  this  definition  includes  both  productive 
and  nonproductive  episodes . This  is  necessary  because  monkeys  have 
cheek  pouches  and  can  conceal  or  retain  vomitus . Thus  it  is  intended  to 
indicate  when  the  subjects  are  distressed  beyond  the  stage  of  nausea. 
Appendix  A illustrates  this  definition  by  graphically  portraying  the 
temporal  nature  of  the  emetic  responses  for  subjects  in  the  repository. 
Sources  of  variation  evident  in  these  figures  include: 

a.  Whether  or  not  a subject  will  experience  an  emetic  episode. 

b.  Number  of  episodes  experienced  if  there  is  emesis. 
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Kind'.,  of  ep  . nodes  ex:  • : oneed  (productive  vs. 

luo t ive)  . 

. Durnt i on  of  episode:  . 


« . ! ir.e  cl*  occurrence  or*  cuch  episode. 

DISTRIBUTION  St  FC I K1 CATION  (Studies  2,  13,  IQ,  27,  10,  VO 

Methods 


: f.e  approach  used  in  modeling  the  prodromal  syndromi  is  that  used 
it.  the  analysis  of  reliability  and  life  data.  The  reliability  of  each 
snl/et*  is  the  infrequenc.N  with  which  lie  experiences  emesis.  Cur  goal 
was  to  estimate  the  probability  of  a failurt — an  emetic  episode  as  well 
as  the  time  until  it  occurs,  and  the  duration  of  the  episode.  The 
decision  tree  in  Figure  1 illustrates  this  approach.  The  origin 
represents  immersion  in  the  radiation  field.  The  probability  of 


it  takes  t minutes  to  occur,  and 


experiencing  a first  episode  is  p 
will  last  d minutes.  Ihe  probability  of  a second  episode  is  p ; it 

minutes  following  the  first  episode,  and  will  last  dominates, 

In  this 


occurs 


L’his  process  continues  until  the  last,  episode  when  p,  . -0 

. (Asd  ' 


manner,  the  prodromal  syndrome 


For  a given 
study  anc  let  x| 
An  estimate  o 


.•  t 


udy , let  n equal 
be  t he 
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viewed  as  a stochastic  process. 

total  number  of  animals  in  the 


; he 


number  of  emetic  responders  during  episode  i, 
p.  = x /n.  Confidence  intervals  for  p can  be 
onstrueted  from  the  binomial  distribution  (of.  ref.  If).  1 


As  stated,  the 


.t. 


c.  e 


fin  it  ion  oi 
upon,  the  interpretation  he  wants  t 
o'  animals  that  have  exactly  i emt 
Feint 


l s 
nke 


vague . 
Let 


'.'tie  modeler  must  decide 


y represent  t tie  number 
ic  responses . y x - x 
parameters  tire  suggested  by  the  additional  estimates: 


P(  - y,  ::  (x.  -xi+t)/n;  O 
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. :.e  approach  to  dot  ora  ini  r*  t,  and  d(  war.  to  consider  then  as 
rat  ion  variables  T and  2,  via  tie  distributions  need  specification.  Wit! 
sna.ll  sample  sites,  custonary  X‘  -goodness-of-fit  tests  are  inappro- 
! " i ; * t for  di  st  r i but  lor.  dot  emi  nation  as  are  the  Kolmogorcv-Sni  rr.ov 

torts  which  require  n priori  knowledge  of  parameters.  We  examine. i 
two  of  the  most  frequently  occurring  distributions  in  life  testing--t he 
exponent  ial  ar.d  the  Wei  bull.  These  densities  are  given  by: 

i 1 \ 


for  the  exponential,  and 


6 \ >■  I 


o X - p ; p.  b ■ Q . P 0 

for  the  Weibiill . With  pi  = 1 and  p = 0,  one  can  see  that  the  exponential 
is  a special  case  of  the  Weibull.  To  test  the  adequacy  of  the 
exponential , we  used  a modified  K2  technique  given  by  Lillieforr,  (5)  for 
samples  3 and  by  Finklestein  and  Schafer  (3)  for  samples  = 2 . A test 
derived  by  Mann,  Ccheuer,  and  Fertig  (9)  was  used  to  test  the  adequacy 
of  a two -parameter  Weibull  distribution  ( p = 0)  when  the  exponential 
distribution  failed.  All  testing  was  done  at  the  ■<  = .01  level. 
Following  selection  of  a distribution,  maximum  likelihood  estimates  were 
used  to  estimate  exponential  paramet  ers  •,  Mann's  (7,  P)  best  linear 
invariant  estimates  were  used  for  I'  and  6 estimates  in  Weibull 
distribut ions . 

Results 

Tables  2 ar.d  3 classify,  respectively,  by  study  and  episode  the 
proportion  of  emetic  responders  and  the  9‘*  confidence  intervals  for 
these  proportions.  The  width  of  the  confidence  intervals  in  Table  3 i.‘ 
a function  of  the  sample  size  and  the  •>  -level  considered.  Tesp'iie 
their  width,  meaningful  data  can  still  be  gleaned  by  examining  the  end 
points  of  these  intervals.  For  instance,  the  lower  limits  for  p in 
studies  2 and  27  indicate  that  at  least  1*77  and  1*3?  of  those 
respectively  exposed  to  such  done  profiles  (and  associated  experimental 
parameters)  would  experience  the  prodromal  syndrome.  Information  of 
this  nature  could  thus  be  used  to  make  best-case  inferer.oes.  By  the 
same  token  upper  limits  for  in  study  30  indicate  that,  at  most,  hit 
of  those  exposed  to  these  experimental  conditions  would  experience  3 
episodes,  while  upper  limits  for  p -p^,  indicate  that  a maximum  of  23? 
will  have  U or  more  emetic  episodes.  ' Worst-case  estimates  can  thus  be 
made  from  these  upper  limits.  Both  best-  and  worst-case  estimates  from 
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in  these  estimate.- , "he 
1:  creasing  the  ■ > -level  . 


init. 


i r.e  wan'  e . 
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. ■ Ldenci  . nt  erv  ■ • uld  be  shorten* 


Tab  ’ es  •.  and  c- 
.••.■'.voted  to  describ 


summar*  ze  by  study  ana  episode  the  distr ibut  ions 
e T and 


( based  upon  ■>  = .01.  In  all  ( IT  of  IT' 
wnere  the  distributions  could  be  identified,  the  distribution: 

*‘‘'.e  • failure  (T-  ) were  exponentially  iistributed.  This  distribution 
suggests  that  when,  th.e  prodromal  syndrome  occurs,  one  car.  char-ict e-- 1 ~e 
when  it  will  happen  as  a function  of  A which  is  dependent  u: on 
experimental  conditions  like  hours  fasted,  radiation  type,  and  th.e 
neutron-lamina  ratio.  Estimates  of  A are  informative,  end  they 
represent  the  mean  tine  between  failures.  In  all  cases,  we  find  that  it 
takes  longer  for  the  first  episode  to  occur  than  for  succe  sive 
epi sodes . 


In  13  of  IT  cases,  the  exponential  distribution  was  an  adequate 
descriptor  of  episode  duration  tines  at  the  .01  level.  The  Weibull  was 
able  to  describe  an  additional  0 cases.  The  remaining  2 cases  were 
m :thtr  Weibull  nor  exponential  ar.d  are  unknown.  With  p = 0,  we  can 
describe  lc  of  IT  cases  of  D, 's  distribution  as  belonging  to  a two- 

parameter  Weibull  family.  This  characterization  suggests  that  episode 

duration  tines  can  be  described  as  functions  of  p and  J>  which  depend 
upon  experimental  conditions.  Thus,  the  sane  mechanism,  appears  tc  bi  at 
work  in  each  study  and  episode,  but  differing  only  in  intensity. 
Estimates  of  A (for  the  exponential  subclass)  are  near,  duration  times 
for  each  episode.  Study  2 with  the  greatest  incidence,  the  most 

episodes,  ar.d  th.e  shortest  tine  to  the  first  episode  appears  to  have  the 

greatest  duration  times  or.  the  average. 


Comment s 

The  diversity  of  experimental  parameters  suggests  a number  of 
otic  es  K be  tested.  For  a given  dose,  they  include  the  hypothesis 


a.  lie  difference 
vs.  a double  pulse,  or 

V.  2 difference 

c.  lie  difference 

d.  No  difference 
Naive) . 

e.  lie  difference 


between,  the  same  dose  delivered  as  a single  pulse 
a single  pulse  vs.  n continuous  dose  rate,  etc. 

bet  wee:,  different  dose  rates. 

between  fasting  animals  IT  h,  t h,  or  1 h. 

between  behavioral  tasks  (FIT,  D+ , D-,  Landolt , 

between  .2^,  .q0,  or  .7?  neutron-gamma  ratios. 


With  hindsight , one  would  like  to  test  each  of  these  hypotheses. 
Giver,  the  current  situation,  none  car.  be  examined  because  no  two  of  the 
several  studies  differ  in  only  a single  independent  variable. 


0 


an  tixe  to  onset 


Ain't  tier  I’oiis  i derat  ion  ec.Td  he  to  shorten  t tie  range  of  total 
!!■  i depigast  r to  doses  t'or  a given  study.  In  study  K! , li  animals  in  tin* 

1 000-rnd  do - e rmu’c  were  deleted  from  cons Iderat  ion  in  distribution 
spec i fi cat  1 ons . This  l . • :'t  7 nnimals  in  the  dose  range  of  CliUO-.'dOC 
rads . It  is  preferaMe  to  make  I nferenees  over  narrow  dose  ranges,  say 
>0  rads,  rut  her  than  over  a 550-rad  range.  This  variability  in  dose 
rany.es  affects  the  precision  of  p.  , p.  , p , \ , f>  , and  |' . 

Additionally,  these  estimates  can  be  improved  by  increasing  tin'  sample 
sices  in  each  study. 


Kl'  PKTKKM1  NATION:!  (Studies  10,  h 1 ) 

Met.  hods 

The  '.T^  number  is  the  effective  dose  for  which  r'0'  of  a population 
will  have  an  emetic  episode.  A literature  seareh  showed  that  estimates 
exist  for  humans,  but  none  could  be  found  for  monkeys.  We  therefore 
sought  to  fill  this  void  t »>  aid  In  monkoy-t . -man  extrapolations.  The 
up-and-down  method  (l,  .')  was  selected  over  probit  techniques  because 
accurate  estimates  could  be  obtained  with  smaller  samples.  The 
procedure  was  to  administer  doses  at  equally  spaced  increments.  If  the 
primate  vomited,  the  increment  was  lowered  one  unit.  if  the  primate  did 
not  vomit,  t lie  dose  was  raised  one  unit  . The  increment  selected  was  r'0 
rads.  Assumptions  upon  which  the  procedure  is  based  include: 

a.  Pat  a is  normally  distributed  with  common  variance. 

b.  i’ub.lects  that  do  not  vomit  would  have  vomited  at  a sufficiently 
higher  dose. 

Advantages  of  this  procedure  are  that  it  concentrates  testing,  near 
the  mean  and  can  save  in  tilt'  number  of  subjects  tested. 

Kesul t s 

The  il\(.  for  unrestrained  man  is  reported  to  be  ll'l  rads 
(log-normal  distribution!  or  . ' 1 i»  rads  (normal  distribution)  (b). 

1st  imated  Kl\(,  for  naive,  stat  ionary,  chair-restrained  rhesus  monkeys 
was  kilt'  rads.  Naive  monkeys  that  were  oscillated  forward  and  backward 
*_  ''  t o IP  (pitch  axis  on  tilt'  ITU’)  from  horizontal  had  an  KPV  of  dt'P 
rads,  rose  rate  was  .’0  ra.ir  min  to  simulate  common  radiation  therapy 
dose  rates  from  which  human  Kl^est  l mates  were  derived.  Kxper iment al 
results  art'  given  in  Tables  ('  and  7.  Large  and  small  sample  Kl^, 
estimates  are  given  in  Tables  0 and  0 . Note  the  close  agreement  of 
estimates  by  both  large  and  small  sample  techniques.  A comparison 
between  stationary  and  motion  groups  Is  in  order.  This  is  accomplished 
via  a large  sample  z-  test  . 
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TARI.K  6.  STATIONARY  TRKATMKNT  TATA 
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TABLE  7 . MOTION  TMKATMVNT  DATA 


Doe# 
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000 


0* 
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0 
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TABLE  6. 

LAME  SAMPLE  KSTIKATI 

■5  (cf.  ref. 

X) 

Group  N 

H A 

®S0L 

( rada) 

(raila) 

0 

(rada ) 

Stationary  8 

T 50 

U U6 . U 3 

27.01* 

Motion  3 

3 50 

258.33 

I8.6I1 

TABLE  9.  SMALL  SAMPLE  ESTIMATES  (cf.  ref.  2) 

Group  H N * d KD30 

(rads)  (rede) 

Stationary  15  l1*  50  ^5.63 

Motion  66  50  258.1*5 


Z = 


v 


x — x 

stations  y motion 

s2  s2 

stationary  motion 


N 


station  aiy 


N 


motion 


44b  43  - 258  33 


5 514 


2 

L 27  04  ^ +^18.64) 


and  is  significant  at  the  5%  level.  We  conclude  that  the  presence  of 
metier  can  trigger  vomiting  nt.  lover  doses. 


Comment s 

'the  up-and-down  technique  uses  substantially  fewer  monkeys  than 
probit  techniques,  arid  should  be  considered  in  radiation  studies 
designed  to  study  effects  of  fasting  on  El\(1  and/or  the  combined 
effects  of  radiation  plus  stressors  other-  than  motion. 

Up-and-down  step  siren  should  be  about  one  standard  deviation  in 
magnitude,  which  approximates  70  rads  for  naive,  stationary  monkeys,  and 
*0  rads  for  naive  monkeys  experiencing  motion  stress. 

This  monkey  IlD^  emesis  data  indicates  a difference  in 
radiosensitivity  between  man  and  monkey  (EI\p  - l8l  rads  for  man  and  M6 
rads  for  partially  restrained  monkeys).  Monkey-man  extrapolation  is 
apparently  not  1 to  1,  and  deserves  further  study. 


II 


DOFF  RATE?  (Studies  VI.  1,  11.0) 

Met  hods 

study  U involves  subjects  receiving  approx  '.mutely  the  same  total 
midepignstr ic  dorr  at  1 of  0 different  dost  rate:.  Mid  1 of  ? different 
fasting  schedules.  1’her.e  fasting  schedul  or  were  designed  to  tent  the 
null  Ivypothoses  that 

(a)  Fasting  does  net  affect  emesis, 

tb)  Pore  rates  do  not  a!  feet  ernes  i a . 


Ml  depignst  r i e doser  ranged  from  M7  t >•  ' 't  rads.  "able  Id  summarirt  a 
the  data  collect  ed  ton;  each  of  0 separat  variable  . Onset  !•  tin  t . Ime 
in  minutes  to  the  first  episode,  with  time  aeeum.u at  i ng  from  the  moment 
the  x-ray  machine  vns  turned  ot  . Offset  tr  the  time  in  minuter,  at  t he 
conclusion  of  the  last  episode.  Total  duration  is  the  total  time  in 
minutes  that  the  sub, loot  experienced  both  productive  and  nonpreduct i ve 
emetic  episodes.  The  period  between  onset  and  offset  times  is  defined 
ns  the  sub.lert 's  vulnerability  period.  The  definitions  of  the  remaining 
variables  fire  self  explanatory.  Figure  0 allows  mean  plots  for  this 
data . 

Resu 1 1 a 

Hypothesis  testing  was  accomplished  by  a two  factor  analysis  of 
variance.  In  all  cases  there  were  no  significant  dose  rate  x farting 
tnt eract ions . significant  dose  rate  effects  (p.  .1)  were  Identified  foi 
onset  ami  offset  times,  duration  time  ot  productive  episodes,  ns  well  as 
the  number  of  productive  and  nonproduct lve  episodes.  Significant 
fasting  effect a (p-.l}  were  found  for  on net  times  and  nonproduct i ve 
duration  times.  Tukey's  (of.  ll)  multiple  comparison  procedure  detected 
the  following  differences  at  the  10*1  level  of  ;!  gni  f icar.cc . 

a.  Onset  and  offset  times  for  0 rads  min  .lose  rate  were  longer 
than  both  the  ‘>0  rad  a /min  and  17*'  rads  min  rates. 

b.  Offset  times  at  1 £ rada/mln  were  logger  than  at  IT1'  rads  min. 

c.  IVoduct  lve  duration  tines  for  1.'  rads/min  were  longer  than  at 
rads/min. 

d.  The  number  of  productive  episode  a*  1 ' rads/min  wore  greater 
( luui  at  r*'  rads, 'min. 

e.  The  number  of  nonproduct  l ve  episodes  at  rads  ’min  were 
greater  than  at  1.'  rads/min. 


Number 

ncnproduct  1 ve 
cpl  no.le* 


, K RATE  DATA 


Oose  r»It  l rada/mln) 


103.1*2 
101. 5* 


127.1*6 
1^1*  .93 


U .3 
63.38 


6.68 

12.8. 


6.68 

7.1*7 


32.70 

52.25 


51.6? 

31.5S 


110.6!* 

82.55 


90.80 

35.32 


86.01. 

39.10 


30.3 

3.7T 


0.93 

3.1!* 


8.3s 

3.77 


8.25 

3.11* 


8.35 

0.0 


0.68 

0.0 


107.00 

1*3.13 


107.08 

121.60 


73.1*7 

3.38 


0.0 

7.00 


V « ;■  f. 


Hour* 
f SSted 


Onset 


06.17 

136.67 


Offset 


122.00 

11*0.97 


Vul nerabllity 


35.91 

?1*.0I* 


Total 

duration 


2.87 

1*.30 


Productive 

duration 


2.37 

1*.30 


N-n  product  Ive 
durnt 1 on 


number 

episodes 


Sunder  productive 
epl 3 odea 


141*  .07 
97.17 


121.62 

11*3.05 


21*.  16 
71.01* 


10.29 

3.1*0 


11.65 

5.13 


0 . 36 

3.i*0 


92.72 

T2.07 


81.63 
5l*  .23 


03.15 

121.15 


121.76 

66.05 


0.1*3 

1*0.13 


1*0.03 

12.72 


0.1*1 

0.82 


I*  .78 
3.95 


0.60 

0.0 


0.1*3 

0.82 


50.50 

67.55 


168.21 
10*..  37 


11.65 
s .13 
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Note  the  lack  of  a clear-cut  pattern  in  significant  findings.  For 
. nstunee,  if  onset  and  offset  times  for  6 rads/min  were  significantly 
1 nger  than  both  the  r0  rads/min  and  175  rads/min  rates,  one  would  have 
ilsc  expected  them  to  be  longer  than  the  92  rads/min  rate.  Similarly, 
if  the  number  of  nonproductive  episodes  at  25  rads/min  were  greater  tlini. 
at  12  rads/min,  then  one  would  have  expected  the  number  of  nonproductive 
episodes  at  50  rads/min,  °2  rads/min,  and  175  rads/min  to  be  greater 
than  the  number  of  nonproductive  episodes  at  the  17  rads/min  rate. 

These  inconsistencies  can  be  attributed  to  the  2 replicates  in  each  cell 
of  the  data  matrices  (see  Table  10).  Too  many  questions  were  asked  with 
too  little  data. 

Comments 


1 

I 1 


Onset  and  offset  times  can  be  measured  from  three  different 
starting  points:  when  the  machine  is  turned  on,  when  it  is  txirned  off, 
or  when  tome  specific  dose  is  delivered  regardless  of  when  the  machine 
was  turned  on  or  off.  When  onset  is  measured  from  machine-on,  a 
mechanical  bias  is  introduced  since  lower  dose  rates  take  longer  to 
deliver  an  effective  dose  of  radiation  to  the  subject.  In  the  above 
data,  if  the  time  to  deliver  l‘!x6  rads  (FD,0  for  emesis)  is  subtracted 
f : om  all  onset  and  offset  times,  then  onset /offset  times  would  measure 
the  lag  time  between  the  delivery  of  an  emetic  dose  and  the  time  of 
appearance/disappearance  of  emesis.  If  dose  rate  has  a physiological 
effect  on  onset/offset  times,  then  those  lag  times  should  change  at 
various  dose  rates. 

Subtraction  of  the  time  to  deliver  Mo  rads  eliminated  all  evidence 
of  dose  rate  effects  on  onset  and  offset  times.  For  example,  onset 
times  measured  from  machine-on  were  107.0  min  at  6 rads/min  and  1x8.7  min 
at  P’r'  rads/min.  These  are  statistically  significant  differences. 

Onset  times  measured  from  delivery  of  1*1x6  rads  were  10.6  mil',  at  6 
rads  min  and  1*6.2  min  at  175  rads/min  and  are  not  statistically 
different.  This  implies  that  the  differences  detected  in  study  31  were 
mechanical  ones  due  to  dose  rate,  and  that  the  phyr i ologicnl ly  important 
event  (between  6 rads/min  and  17r'  rads/min)  was  not  dose  rate  but  the 
accumulation  of  an  effective  emetic  dose. 

Next  , we  use  Type  1 ( •>  ) and  Type  li  ( •'  ) errors  to  discuss  the 
implications  of  testing  t dose  rates  ai  d .'  fasting  schedules  at  once. 
Recall  that  a i'ype  1 error  is  the  probability  of  rejecting  the  null 


hypothesis,  when  it  is.  true.  For  this  data  set 


Type  T error  in 


committed  when  wo  find  either  significant  dose  rate  or  fasting  effects 
when  there  should  be  none.  A Type  IT  error  is  committed  when  one 
accepts  the  null  hypothesis  when  it  is  false.  For  this  data  set,  a I'ype 
II  error  is  committed  when  one  finds  neither  significant  dose  rate  nor 
fasting  effects  when  there  should  have  been  some . We  believe  this  may 
be  the  case  in  the  present  situation. 


Ar.  nlternat  ;ve  to  measuring  Type  I T error is  t . examine  the  power 
of  the  test  where  power  is  l-  if  . In  hypot  her.  5 r.  testing  the  investigator 
specifies  the  .1  -level  he  can  tolerate  an<!  t hen  chooses  the  sample  si 
to  maximize  the  power  of  his  test.  In  this  way  or.,  minimizes  the 
probability  of  committing  a Type  II  error.  Figur.  <-r- , respectively, 
show  power  curves  for  testing: 

n.  Both  the  hypothesis  of  no  dose  rate  « f foe*,  among  b dose  rate: 
and  the  hypothesis  of  ro  fasting  effect  among  far*  ir.g  schedule:  . 

b.  The  hypothesis  of  no  fasting  effect  •>•:.  eg  fasting  schedules. 

o.  "'he  hypot  ties i s of  no  dose  rate  effect  ar  er.g  t dose  rates. 

Delta  (A)  is  the  maximum  d • • ference  be twee:  dose  rates  or  fasting 
schedules  that  the  invest  igator  r..ij  wish  to  detect  . >igma  ("  ) is  the 
estimated  standard  deviation  for  each  variable . Dr.ch  estimate:-  (the 
^ MSF  from  the  ANOVA's)  are  summar i zed  in  Table  . 

To  illustrate  the  use  of  Figures  !-*• , suppose  ar.  investigator 
wishes  t let  cot  10-m.ir.  onset  time  iifferenotr.  in  dose  rater  and  fast  ir.g 
schedules  at  the  s?  le\el  of  rignlf : ear.ee.  In  such  a case  delta  sigma  = 
10  .ID. do  or  1.1.  From  Figure  o we  see  that  the  probability  of  a correct 
decision  (the  power)  is  iOf  for  the  experiment  ns  it  wns  implemented 
with  Di)  subjects.  By  increasing  th.c  sample  size  by  a factor  of  five  t, 
IDO  subjects,  a more  realistic  power  level  of  hUf  results . If  the 
invest  igator  were  willing  to  omit  the  question,  of  o.ose  rate  effects  and 
concentrate  only  on  the  question  of  fasting  effects  (with  •>  r .0rO,  the 
same  .'k  animals  would  yield  n power  level  in  excess  of  QO*  (cf.  Fig.  •<  1 . 
On  the  other  hand,  if  one  were  willing  to  use  to  animals  to  test  the 
Ivypothesis  ( >»  = . 0C' ) of  effects  among  the  6 dose  rates  (while 
sacrificing  an  answer  to  the  question  of  fasting  effects',  the 
probability  of  making  a correct  decision  would  be  ’lit  (of.  Fig.  ^ ' . 

Figures  and  the  above  show  that  there  are  a number  of  ways  to 
increase  the  probability  of  making  a correct  decision.  They  include: 
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Figure  4.  Minimal  (worst  case)  power  curves  for  testing  both  t lie 

hyjiothesis  of  no  dose  rate  effect,  among  6 dose  rates  and  the 
hypothesis  of  no  fasting  effect  among  2 fasting  schedules. 
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Figure  !*.  Pover  curves  for  testing  the  hypothesis  of  no  fasting  effect 
among  2 fasting  schedules. 
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gure  r- . Power  curves  for  testing  the  hypothesis  of  no  close  rate 
effect  among  6 dose  rates. 


Accurate  estimates  of  •»  are  essential  to  ensure  that  operationally 
significant  differences  can  he  detected.  In  the  above  example  30  min 
represented  a difference  of  1.3  standard  deviations  in  onset  times.  In 
real -world  applications  it  might  turn  out  that  smaller  differences  in 
onset  times  (say  10  min)  are  relevant  differences  for  detection.  This 

would  entail  even  larger  sample  s i ses . We  recommend  at least  10 

subjects  treatment  group  for  intergroup  comparisons. 

MIXED  RATES  (Studies  3b,  37,  38,  Uo) 

Met  hods 

Studies  !b,  37,  3t',  and  40  represent  a departure  from  previous 
experiments  in  that  they  call  for  scenario-specific  profiles.  Table  10 
summarizes  their  experimental  parameters.  All  four  studies  in  this 
section  use  mixed  dose  rates  protracted  over  relatively  long  exposure 
periods  (0-T0  h) . These  studies  must  be  viewed  as  pilot  studies  due  to 
their  small  sample  sizes  and  their  unique  dose  profiles. 


01 


TABLE  12.  MIXED  RATE  EXPERIMENTAL  PARAMETERS 


Sequence 


Dose  rate 


Durnt Ion 


Tir.e  dose  rate  Is  given 


Study  36  (Task  « D-,  work.  period  * 11/12  h)  ar.d  Study  37  (Task  = D+;  work  period  = 11/12  h) 


27.6  rads /min 
97.8  rads/h 
9.9  rada/h 
27.6  rads/nin 


2 min 
.62  h 

10  h 

3 min 


Suanary : 


No . dose  rates 

U 


Total  exposure  Derlod 
12  h 


Total  midepigastrlc  dose 
302  rads 


Study  36  (Task  = PEP;  work  period  = 9/72  h) 


1 

13.5  rads/mln 

10  mic 

0 h 

2 

10  rads/h 

10  b 

10  min 

3 

1.05  rado/b 

62  b 

10  h 

Summary : 

No.  dose  rates  Total 

3 

exposure  period 

72  h 

Total  r.idepi gastric  dose 
300  rads 

Study  39  (Task  = PEP  ♦ MART;  work  period  « 10/10  h) 


1 

50  rads /min 

2 min 

0 h 

2 

17.5  rsds/min 

2 min 

1 h 

3 

7 rads/mln 

15  min 

2 h 

U 

1000  rads/mln 

2 min 

T .5  h 

Summary : 

No.  dose  rates  Total 

U 

exposure  period 

10  h 

Total  midepigastric  dose 
lLl*0  rads 

22 


ftudies  36  and  37  used  the  3ane  profile  to  study  the  effects  of 
both  positive  (D+)  and  negative  (T-)  reinforcement.  Animals  vere 
trained  to  perform  identical  tasks  vith  different  rewards  and  vere 
expend  to  identical  radiation  profiles.  One  group  received  food 
pellets  for  a correct  response.  The  other  group  received  no  food  but  a 
shock  vher.  they  made  an  incorrect  response — negative  reinforcement. 

Table  13  shows  the  number  of  animals  in  each  study  that  had  at  least  one 
emetic  episode. 


TATLE  13.  »2*IIC  RESFONSES  IN  STUDIES  36  AND  3? 


Shock 
Study  36 


rood 
Study  37 
IH 


No  ciesla 


No.  of  subjects 


Results 

Fisher's  exact  test  (cf.  10)  found  no  difference  in  the  proportion 
of  emetic  responders  between  the  two  studies  at  the  lOt  level  of 
significance.  This  finding  is  a result  of  the  sample  size  selected.  If 
r.c  shock  animals  had  vomited,  Fisher's  exact  test  would  have  found  a 
significant  difference  between  studies  at  the  .05  level.  Attention  must 
be  paid  to  the  questions  to  be  answered  and  the  sample  size  selected  to 
answer  these  questions . 

It  would  be  nice  to  compare  study  3?  with  studies  36  and  37  as 
these  studies  involved  30C-rnd  total  midepigastric  doses.  However,  as 
enumerated  in  Table  lb,  a number  of  differences  between  these  studies 
are  based  solely  on  selected  experimental  parameters. 


TABLE  1U.  EXPERIMENTAL  PARAMETERS  FOR  STUDIES  3b,  37,  AND  38 


Studies  ?6  ar.d  17 


Discrete  avoidance  (lever  pressing' 

Stationary  subjects 

different  dose  rales,  durations, 
and  times  of  administration 
U work  sessions  vlth  TO-sdnute 
rest  periods  between  sessions 


1 different  dose  rates,  durations, 
and  tlr.es  of  administration 
4 1 'T-h  work  sessions  vith  rest 
periods  of  1/t  h or  It  h between 

sessions 


I 


f 


\ 


L 


Even  if  the  sample  nlaes  wore  nut' fir  lent  l.y  large  to  permit  valid 
numerical  comparlHons  between  those  studies,  one*  would  be*  hnrd  pressed 
t o i! t.ate  the  cause  of  any  significant  differences  det.ee t e*d ; should  t hoy 
be*  attributed  to  tank  differences,  motion  di  f 1 rrrnces , dose-rate 
d I ffereneeB,  or  work-rest  sehedul  es  . 


In  study  19,  subjects  were*  trained  to  perform  t,-.  .. 
simultaneously — the  continuous  PER  (pitch  axis  only)  and 
MART.  The*  done  rnten  differe'd  from  the  * ther  ntudlen  in 
(see  Table*  It);  in  addition  the*re»  were*  no  re*st  periods, 
differed  from  other  studies  throughout  (tie*  itq  > i;  i torj  in 
only  P-task  e*x  pe*i*  ime'ii  t that  war  i mpl  eme  nt  e*d  . tnterstudy 
we*re*  not  pe'ssihle*  because  of  the*  variety  of  e*xpor i me  tit  a.l 


citliere*nt  tanks 
the*  discrete 
this  section 
This  study 
that,  it  is.  the 
rompa r i sous 
pariune't  eta  . 


t'e  ’Mime*n  t 


M.xrd  rate*  experiments  may  provieie*  a practical  are -a  of  rise  arcl  : 
a.  'Peed  inn  for  orde*r  effects  In  dose  sequences, 

h.  Testing  foe*  time-  differences  for  the*  same  done  sequence*, 

e.  Tenting  for  time  differences  for  done  rate*  sequence*n. 

d.  Tenting  for  tank  differences  for  the*  name  dose*  profile. 


Each  of  these  arean  can  he  systematically  examined  and  tlms  provides  ei 
foundation  upon  which  general  statements  can  he  made.  The  alternative* 
is  to  continue*  looking  at  scenario-specific  profiles  which  can  only 
answer  wheit.  might  happen  under  a given  set  of  conditions.  Such  an 
approach  provides  little  or  no  information  U field  commanders , if  the 
reality  of  their  situation  doe:  not  Closely  resemble  the  experiment  a.l 
r, conn r ! o report  el . 

DESCRIPTIVE  (Studies  3.1,  3.2,  Pb,  Pfl,  IP.  tit.  3b) 

This  section  contains  stud  lea  from  the  repository  which  could  not 
he  classified  under  any  of  the  other  headings.,  oimilarly,  they  are  also 
pilot  studies. 

Pub  .loots  In  studies  3.1  and  3.2  were  trained  lo  perform  the  handol  t 
ring  visual  acuity  task.  Table  IS  shows  subjects  In  the  two  studies 
were  exposed  to  1 of  doner.,  1 of  3 fasting  schedules,  of  P cent  lulled 
dose  raf.es.  and  the  same  n : ) ratio.  Multiple  experiments  1 parameters 
coupled  with  small  sample  sires  preclude  the  testing  of  hypotheses  of 
differences  between  doses,  dose-  rates,  and  fast  i ng  schedules  for  11  i: 
behavioral  task  and  this  neut ro a -gamma  rat  in. 


P)» 
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;t  mly  l.l 

IVm»  (ISils.*,  n ) * .S) 
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6 
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No.  of  t*o t h Hi  Kh  1 t'l\  lit  0\ \tOi  l'l  s \ 

Wtte  h l « 1 C 1 >* 

No  !‘nr u l « 0 

Total  ' 2 lo 

iti  I't  ndy  all  subjects  were  on  the  same  1'aat  tng  schedule  (if 
dose  rate  (pulse),  neutron-gamma  ratio  (.6r>),  and  tank  (PFF).  They 
be  grouped  into  It  dose  ranges  as  shown  in  Table  It'.  Biologic 
variability  and  amall  sample  sizes  lead  to  a relatively  flat 
dose-response  curve.  Fines  is  did  not  occur  at  the  highest  doses  (.101 
and  .W'lt  rads)  tested. 
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2 
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Study  28  was  similar  to  study  .'6  in  that  tin*  tusk  (PEP),  doae  rate 
(pulse),  neutron-gamma  ratio  (.'>',  and  fasting  schedule  (l  h)  were  held 
constant  for  all  subjects.  The  small  sample  Hive  spread  over  many  doses 
prevents  definition  of  a dose-response  curve  (Table  17).  This 
experiment  did  serve  the  purpose  of  dose  ranging , by  identifying  the 
occurrence  of  emesis  in  (tie  Slip-  to  1110-rad  dose  interval. 


tahi.k  17.  wtFTic  KmvNsrs  in  study  eft 


pour  (mill 

UViUr,  n y 

- .5) 

i .■  ■ • 

Site  ASp 

9k  5 | g 

\°  lo-i  no 

Total 

ftu*al  » 
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? 
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6 

No  coueals 

3 

0 

0 

l 

it 

No.  of  subjects 

3 

3 

2 

10 

In  study  32  all  naive  animals  were  subject  to  the  same  experimental 
parameters:  fasting  = 1 h;  n:>*  l /U ; dose  rate  pulse;  dose  range  « QSP 
- 1130  rads.  This  study  Is  uni  .pie  in  that  only  of  12  subjects  vomited 
in  thin  relatively  high  dose  range  within  the  h postexposure  period. 

Study  ilt  compared  2 doses  t ‘>00  and  1000  rads'  using  naive  animals 
exposed  to  a 5/1  neutron/ gamma  ratio.  As  shown  in  Table  li', .the  small 
sample  siv.es  and  different  fasting  schedules  precluded  the  making,  of  an,\ 
rigid  statements  except  that  emesis  can  occur  at  either  ''00  or  1 0CK 
rads . 
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10 

Total 
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ta.'k ) 
lose 
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sly  <'  was  allot  he l ranging  experiment  involving,  t .1  1 ffeeent 

who:  e ubjeot:.  were  oxpe  od  to  the  same  experimental  pnraiuetera  ta 
•mi  ik  ratio;  n 1 h faat.ing  schedule,  and  the  ITT  plteh  and  roll 
''able  19  annunar  1 the  , ‘ourrenee  of  emesis  at  the  various 
Onset  time!  were  earlier,  and  the  average  number  of  episodes 
rent er  for  the  lOCV-rnd  dose. 
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The  data  from  10  studies  i t,  K' , Id,  .'7,  .’0,  U,  11,  ts,  19) 
were  arbitrarily  grouped  together  to  provide  a very  gross  doner t pt ion  of 
the  emet le  response  following  radiation  (Table  TO).  The  group  total  was 
>'h)  monk.ys  representing  all  aspeots  of  experimental  variability  suoh  as 
different  doses,  dose  rates,  radiation  quality  (x,'.  n:>),  fasting 
times,  and  tasks.  The  oMeottve  Is  simply  to  provide  future 
experimenters  with  rougl  statist  leal  estimates  to  plan  resout 
oomm  I intents,  hopefully  Ineludtng  fewer  cpiest  Ions  , more  monkey;'  per 
treatment  group,  and  longer  observation  periods. 
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00M1VTKH  PROGRAM;; 


This  section  briefly  outlines  the  computer  software  used  in  the 


analysis  of  this  re 

tory . 

All 

programming  was  done  in  API, 

. Source 

listings  of  appro pri 

a i e code 

are  $ 

;iven  in  Appendix  R.  This  section  is 

not  an  Ar'T  tutorial  . 

Tt  s purpose 

is  to  indicate  by  name  and 

fund  ton 

appropriat e soft  ware 

tools  sc 

' that 

the  repository  can  be  maintained. 

Out  a 

Sets 

Data  set  KMES Tit 

.0770.1, '3 

is  organized  by  study  as  follows 

Description  of  data 

F i e 1 d 

Description  of  data 

Field 

Study  number 

!_s 

Animal  type 

li  1-1*8 

Month 

7-9 

Number  of  males 

50-c>‘ 

Year 

11  -12 

Number  of  females 

s7-tl 

Experimental  site 

1 t-lf' 

Restraint  used  on  subject 

6 i-t'7 

Type  radiation  used 

70-214 

Emesis  data  collected 

t>o_7b 

Neut  roti /gamma  rat  i o 

e~.'d 

Low  dose  (M  i depigast.fi  c ^ 

77-80 

Task 

U - ill 

High  dose  (Mi  dept gnst.r i c ) 

8 3 -8p 

Hours  observed 

1b-  17 

Mode  of  dose  administered 

88-01 

Number  of  animals 

Uo-Ui 

Tlie  function  REPOSITORY  accesses  this  data  set.  It  is  user 
oriented  and  permits  the  novice  to  examine  the  experimental  parameters 
in  ail  studies  currently  cataloged  in  the  computer  repository.  This 
function  can  produce  tables  similar  to  .able  .1  as  veil  as  subtables  of 
any  of  the  experimental  parameters  selected.  For  example,  one  can 
obtain  a listing  of  all  ITT  experiments  or  all  experiments  conducted  at 
White  Sands  Missile  Range. 

The  data  set  ALI. .EMESIS .0770420 , organized  by  subject  and  study, 
contain:  additional  experimental  parameters  on  each  animal  such  as 
fasting  schedule  and  dose  rate  as  veil  ns  what  happened  on  exposure  day: 
the  number  of  emetic  episodes,  onset  and  duration  times,  and  the 
c 1 ass  i f i cat  i on  of  episodes  as  productive  or  v.onproduct  i ve . This  data  is. 
structured  as-  follows- 


Description  of  data 

Field 

Description  of  ,4ut a 

FAshi 

Subject  ID 

1-8 

!\tl  se  s i y. t (min) 

Dl-Db 

Study  No. 

7-10 

Dose  rate  (rads  min) 

28-  31 

Dose  (Midepigast vie ) 

12-15 

Number  of  episodes 

3.3-:3l* 

Hours  fasted 

1 7 -DO 

Onset  t iru  (xxx.xxmin' 

3D -ill 

Radiation  mode 

tCC. 

Duration  time 

li  ?-l,7 

(0  * steady  or  mixed 

rates  1 

Kind  of  emesis 

to 

(l  * single  pulse' 

( 1 • product  1 ve , D * 

ronprod ' 

( 2 * double  pulse' 

fit',  of  contractions 

51-5  l 

d8 


— — •— — 


X. 


<•  pattern  starting  in  column  <6  for  onset  time;.,  duration  time:,  kinds 
ernes  1 s , and  number  of  contractions  repeats  itself  (starting  in  column 
a:  r : ny  times  as  is  necessary  tint  i 1 all  episodes  have  been 


riot  Koutines 


The  graphs  In  Appendix  A were  created  by  the  function  FLOTEMK. 

This  function  obtains  its  data  from  the  ALL .EMESTF .D770U20  data  aet . 
When  called,  this  function  will  ask  the  user  to  define  the  study  he 
wishes  to  draw.  The  function  will  automatically  draw  dose  profiles  for 
cont inuous-  and  single-pulse  experiments.  The  user  must  write  his  own 
software  to  draw  more  complicated  dose  profiles.  PRA33,  PRA3o,  PRA  IT , 
and  DRA30  are  such  examples  for  studies  33,  3b,  37,  and  38, 
respect i vel y . 


The  power  and  the  dose  rate  curves  were  drawn  by  the  functions 
PKTl'rVOP  following  correct  definition  of  X and  Y variables.  These  APT 
functions  are  general  in  nature  and  permit  the  user  to  draw  as  many 
coordinate  axes  with  ns  many  variables  as  one  specifies. 

A graphics  PIRIVTOKI  has  been  created  which  stores  all  the 
illustrations  created  in  this  technical  report.  To  display  any  Item  in 
the  directory  one  merely  needs  to  type  the  word  PIP FI  .AY  followed  by  the 
name  of  the  file  to  be  displayed  in  quotes.  Tills  system  will  save 
dollars  and  time  In  not  having  to  recreate  these  1 1 lust  rat  ions  for  a 
second  time.  A second  advantage  of  this  system  is  that  it  allows 
graphic  overlay's  so  that  different  studies  can  be  examined 
simul tnneoua ly . 

Analysis  Programs 

Analysis  programs  developed  in  this  report  Include  I.ILLYKX  X, 
VKUU’Ll  X,  AM OVA  X,  TFKFP,  and  DFTAT  X.  These  routines  are  user 
oriented  and  will  prompt  the  user  for  necessary  information.  LII.1YF.X  X 
urd  WKiRPld  X,  respectively,  test  tv'  see  it’  data  X has  exponential  or 
Wei  bull  distributions.  AX OVA  X Is  used  to  perform  a two-factor  analysis 
of  variance  on  a three-dimensional  data  matrix  X where  the  first 
dimension  represents,  replicates,  levels  of  the  1'  factor  are  the  second 
dimension,  an.1,  the  levels,  of  the  A factor  are  the  third  dimension. 

kKO  will  take  the  output  of  ANOVA  X and  perform  Tukey's  multiple 
comparison  procedures  on  each  factor.  PP'l’AT  X will  compute  means, 
standard  deviations,  and  other  descriptive  statistics  on  the  data  vector 
\. 

Again,  we  repeat  that  the  purpose  of  this  sect  ion  and  Appendix  H is 
not  to  serve  as  a "how  to"  guide.  Their  Inclusion  is  to  indicate  which 
tools  are  currently  available  in  support  of  this  repository.  There  is 
no  substitute  for  computer  terminal  experience  and  a reading  of  the  R'.'V 
Use t 1 c Repository  log  Rook  to  better  learn  how  t o use  these  tools. 
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Data  on  ionising  radiation  enesi:  , collect  el  from  studio;, 
involving  010  rhesus  monkeys,  has  been  examine,  . "her-«  studies 
contained  a wide  variety  of  experiment a]  parnn.e * < v.  prepared  by 
investigators  interested  in  unique  situations.  h<  tudies  vi  n group  d 
a posteriori  on  the  basis  of  the  analysis  pern.  'to.  by  i'.o  data 
ool leot  od  . 


In  several  studies  the  distribution  of  :;j  • onset  times  and 
episode  durations  could  bo  determined.  In  thoso  i nst  suice;. , * he  onset 
time  of  each  episode  had  an  expedient  ia  \ distribution,  and  tlie  duration 
of  each  episode  vns  either  an  exponential  or  Wei  bull  distribution.  The 
presence  of  similar  distributions  across  many  experiments  suggests  that 
the  same  mechanism  is  at  work  in  each  study  and  episode,  differing  only 
in  intensity. 

One  study  found  the  UDS1  fur  radiation  emesis  to  l o U i*t'  rads  (Co*"  , 
00  rads /min),  and  that  motion  reduced  the  IT v to  .“-8  rads.  Another 
study  (Ah'.  VI  ) demonstrated  lOCt  incidence  of  emesis  at  ‘O'O  rads  ( X-ray, 
6 to  171?  rads/mln). 

Although  designed  to  answer  other  research  questions,  these  studies 
are  good  emesis  syndrome  pilot  studies  in  that  they  refine  experimental 
techniques,  permit  the  training  of  personnel  and  the  development  of 
SCFs,  and  provide  experience  for  contingency  plans.  They  provide  the 
opportunity  to  operationally  define  what  is  to  be  measured;  they 
facilitate  the  selection  of  competing  metrics,  and  suggest  hov  long  the 
measurement  process  should  occur. 

The  effects  of  such  factors  as  dose  rate,  fasting  time,  radiation 
quality,  and  performance  tasks  on  the  radiation  emesis  syndrome  are 
still  undetermined.  The  present  studies  provide  estimates  of  <>  for 
planning  fully  developed  emesis  experiments  so  understanding  of  the 
interrelationships  between  i> -level s , number  of  questions  to  be  answered, 
operational ly  relevant  differences  to  be  detected,  and  the  sample  sine 
selected  can  result  In  cost-effective  experiments  with  a high  payoff  in 
valid  findings.  As  a rule  of  thumb,  we  recommend  at  least  10  subjects 
per  treatment  group  in  making  intergroup  comparisons. 
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A1TENDIX  A 


PICTORIAL  SUMMARY  CP  EMETIC  FINDING? 

Appendix  A is  a graphical  summary  of  emetic  findings  for  the 
studies  listed  in  Table  I.  Stars  2 units  high  indicate  productive  epi- 
sodes. Bars  1 unit  high  indicate  nonproductive  episodes.  Width  of 
bars  indicate  episode  dilation.  Flat  lines  show  no  emetic  activity 
between  emetic  episodes.  Absence  of  bars  and  lines  indicate  no  emetic 
signs  during  the  observation  period  given  in  Table  1.  This  appendix  is 
organized  by  study  numbers  presented  in  Table  1. 
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appendix  b 


COMPUTER  SOURCE  CODES 


> ;...enuix  B is  n listing  of  computer  source  codes  from  which  the 
present  repository  analysis  and  maintenance  was  possible.  The  major 
functions  include: 


REPOSITORY 
FLOTEME 
SETUPMOD 
DIRECTORY 
LILLYEX  X 
WE I BULL  X 
ANOVA  X 
TUKE2 
DSTAT  X 

Appropriate  subroutines  are  also  listed.  The  inclusion  of  this  appendix 
is  to  identify  computer  tools  presently  in  support  of  this  repository. 

It  is  not  a "how  to"  guide. 
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